Causality vs. Ward identity in disordered electron systems 
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We address the problem of fulfilling consistency conditions in solutions for disordered noninter- 
acting electrons. We prove that if we assume the existence of the diffusion pole in an electron-hole 
symmetric theory we cannot achieve a solution with a causal self-energy that would fully fit the Ward 
identity. Since the self-energy must be causal, we conclude that the Ward identity is partly violated 
in the diffusive transport regime of disordered electrons. We explain this violation in physical terms 
and discuss its consequences. 
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Quantum charge transport in random media is gen- 
erally a rather complex problem. Unfortunately, even 
the simplest quantum tight-binding model of disordered 
noninteracting electrons does not provide for an exact 
solution. To find a complete solution of the Hamilto- 
nian with random energies amounts to the determina- 
tion of the spectrum of the Hamiltonian together with 
the spatial extension of the eigenfunctions. If the disor- 
der is sufficiently small the eigenfunctions are extended 
waves enabling diffusion with a metallic conductivity. As 
shown by Anderson in his seminal paper on quantum dif- 
fusion, Ref. 1, strong randomness in the scattering poten- 
tial can, however, destroy long-range tails of the extended 
eigenfunctions and cause electron localization. Localized 
eigenfunctions do not contribute to diffusion and drive 
the system insulating. After years of intensive studies, 
full understanding of Anderson localization has not yet 
been achieved and it still remains a big challenge in con- 
densed matter theory. 

Although there are efficient mathematical tools to 
handle finite random matrices, it is more practical to 
use configurationally averaged Green functions to charac- 
terize solutions of random systems. Green functions have 
proved useful in particular in the thermodynamic limit 
where the volume and the number of particles are infinite. 
Exact solutions are mostly elusive and we have to resort 
to approximations. A mean-field theory of the Ander- 
son model of disordered electrons, the coherent-potential 
approximation (CPA),0 contains only extended states. 
Various extensions beyond the CPA are unable to include 
localized eigenstates either. To take into account local- 
ized states, one has to treat and approximate one- and 
two-electron averaged Green functions simultaneously. 
The one-electron resolvent determines the spectrum of 
the solution, while the two-electron resolvent determines 
the spatial extension of the underlying eigenstates. 

To arrive at a consistent solution with one- and two- 
electron averaged Green functions we have to comply 
with a number of restrictive consistency conditions re- 
flecting analyticity and causality on the one side and 
conservation laws on the other. First of all we have to 



warrant that the eigenenergies of the random Hamilto- 
nian remain real. This is the case if the averaged re- 
solvent G(k, z) = i^s\\_z\ — H~\ 1 |k}\ has no poles 

\ / av 

in the upper and lower half-planes of complex ener- 
gies z. We denote here |k) the eigenfunctions of the non- 
random kinetic energy (Bloch waves), labelled by mo- 
menta k. The analytic behavior of the averaged resol- 
vent is expressed in the Kramers-Kronig relations be- 
tween its real and imaginary parts. If we use the Dyson 
equation to represent the averaged resolvent via a self- 
energy G(k, z) = [z — e(k) — £(k, z)]~ we can write the 
Kramers-Kronig relations for the self-energy 
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where £oo is the constant part of the self-energy (its high- 
energy limit). 

When constructing two-particle functions we have to 
obey besides analyticity also conservation laws expressed 
in Ward identities. Vollhardt and Wolfle proved the fol- 
lowing integral Ward identity in the Anderson model of 
disordered electrons 

£ (k+ , z + ) - £ (k_ , z_ ) = 1 Aft, (z+ , z- ; q) 

k' 

x [G(k' + ,z+)-G(k'_,z_)] , (2) 

where A eh stands for the irreducible vertex in the 
electron-hole scattering channel, k± = k ± q/2, and 
z± = E ± lo ± irj, t] > 0. The original proof of iden- 
tity @ did not reveal the conservation law reflected by 
this identity, but it was shown later that Eq. is a con- 
sequence of probability conservation in the Hilbert space 
of Bloch waves. 

Identity @ plays an important role in physics of dis- 
ordered electrons, since it sets up a condition for the 
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existence of the diffusion pole in the density-density cor- 
relation function. [3, |8j The diffusion pole in the electron- 
hole correlation function is a hallmark for the diffusive 
character of charge transport in disordered systems. 

Neither Eq. (J2J nor Eq. are automatically guaran- 
teed in quantitative treatments of disordered electrons. 
In fact, different conditions must be fulfilled to comply 
with cither Eq. (|TJ or Eq. (J2J. These conditions can be 
contradictory in specific situations and it need not be 
always possible to comply with both equations simulta- 
neously. It is the aim of this Letter to materialize such 
a situation in models of disordered noninteracting elec- 
trons. The principal result of this Letter is the following 
Assertion: If we assume electron-hole symmetry in the 
vertex functions and the existence of a diffusion pole in 
the low-energy limit of the electron-hole correlation func- 
tion having the explicit form 2irnp/[— iuj + Dq 2 }, it is 
impossible to accommodate Eq. © for z + — z_ ^ with 
causality of the self-energy, i. e., Eqs. IjTJl. 

Before we proceed with a proof we first briefly discuss 
the assumptions in this assertion and illustrate in simple 
terms how a conflict between causality and the Ward 
identity emerges. 

To see how causality of the self-energy can be traced 
in the electron-hole irreducible vertex we reduce the com- 
plex frequency difference in Eq. (0) to a purely imaginary 
value, i. e. z + — z_ = 2ir\ and set q = to obtain a special 
case of the Ward identity 
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We have denoted A{£,(z + , z_) = A{&,(z+,z_; 0). 
Causality of the self-energy (negative semi-definiteness 
of its imaginary part) can then be controlled by positive 
semi-definitcness of the self-adjoint part of the electron- 
hole vertex appearing on the r.h.s. of Eq. ©. A causal 
self-energy can be determined from a given positive semi- 
definite electron-hole irreducible vertex. One uses the 
Ward identity to resolve the imaginary part of the self- 
energy from the electron-hole vertex and completes it 
with Eq. (|Ia|l to determine the real part of the self-energy. 

The self-energy calculated from the electron-hole irre- 
ducible vertex via Eqs. © and (|Ia|l is causal and uniquely 
determined. We can warrant fulfillment of the Ward 
identity for imaginary energy differences, i. e., lo = 0, 
in this construction. Behavior of the electron-hole vertex 
for u> ^ is irrelevant for the self-energy. Whether the 
self-energy from Eq. (© satisfies the Ward identity also 
for real-energy differences remains to be checked. 

We know that the mean-field solution (CPA) obeys 
the Ward identity for arbitrary difference of complex en- 
ergies. It is due to locality of the CPA electron-hole 
irreducible vertex. To find out whether Eq. © can 



generally be fulfilled we hence have to analyze nonlocal 
electron-hole irreducible vertices. 

The simplest nonlocal contributions to the electron- 
hole vertex are crossed diagrams. The self-energy is 
causal when the electron-hole irreducible vertex is a (pos- 
itive) sum of products GG* . Perturbation theory up to 
second order for the electron-hole vertex with this prop- 
erty reads 



A eh 
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where the dashed line stands for the local impurity poten- 
tial A = c ((V 2 ) av — (V)a V ) and the oriented solid lines 
are renormalized electron (retarded, upper line) and hole 
(advanced, lower line) averaged propagators. The ap- 
proximate self-energy constructed from the vertex A^ 1 
via Eq. 10 and Eq. (|Ia|l is causal. It does not, however, 
obey the Ward identity © as can easily be checked by 
means of its proof. @ 

To comply with the Vollhardt-Wolfle identity (J2J) up 
to second order we have to supplement the simplest (max- 
imally) crossed diagram with two further diagrams of the 
same order. That is, instead of Eq. © we have to con- 
sider the following electron-hole vertex 
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We can find a formal diagrammatic solution for the self- 
energy corresponding to the vertex A^ 7 : 



(5b) 



The Ward identity holds for T,wi and A^ 7 , but the self- 
energy is not causal. The imaginary part 3£^ 7 becomes 
positive for specific energies and disorder strengths. It 
hence does not fulfil the Kramers-Kronig relation lla|) 
and consequently it cannot be obtained from Eq. I|5a|) via 
Eq. ©. To amend this we have to add an infinite series 
of compensating diagrams (from the so called traveling 
clusters) so that the resulting self-energy representation 
generates a Herglotz function. [Ifj| 

The situation cannot be remedied by adding higher- 
order maximally crossed diagrams, i. e., diagrams with 
scatterings only between the electron and the hole where 
each impurity line crosses all other impurity potentials. 
With each maximally crossed diagram of order n, gener- 
ating a causal contribution to the self-energy, we have to 
add n additional compensating non-causal diagrams. 
The corrected irreducible vertex, compatible with the 
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Ward identity, generates a non-Herglotz self-energy and 
the sum of diagrams is not causal at any order of the 
expansion. 

It is impossible to find a diagrammatic representation 
of a causal self-energy obeying Eq. @ for real frequency 
differences and containing a sum of maximally crossed di- 
agrams. Inconsistency between causality and the Ward 
identity in perturbation theory originates in these nonlo- 
cal diagrams. To prove incompatibility of causality with 
the Ward identity nonperturbatively we need an infinite 
sum of maximally crossed diagrams producing a singu- 
larity in the electron-hole irreducible vertex. This sin- 
gularity, the Cooper pole, is generally an image of the 
diffusion pole in the full electron-hole vertex. The dif- 
fusion pole in the low-energy limit of the electron-hole 
correlation function was proved to have exactly the form 
assumed in the Assertion. Q 

Proof of the Assertion. The diffusion pole in the 
electron-hole correlation function is a residue of a pole 
in a more fundamental full two-particle vertex (Green) 
function Tkk' (z+ , Z- ; q). We denote A eh and A ee the 
irreducible vertices in the electron-hole and the electron- 
electron scattering channels characterized by antiparal- 
lel and parallel propagator lines, respectively. Let fur- 
ther J be a vertex irreducible in both channels. The 
following identity A eh + A ee = I + T then holds. Since 
the vertex / is regular, the pole from the full vertex F 
must appear at least in one of the irreducible vertices. 
Due to the electron-hole symmetry, Af^,(z+, z_; q) = 
A^, (z + , z_; — q — k — k') both irreducible vertices share 
the same singularity. 

The explicit singular low-energy asymptotics of the 
electron-hole irreducible vertex A eh can be obtained from 
the sum of all maximally crossed diagrams. At zero tem- 
perature for q — and k + k'| — ► this sum yields 



kk' ( z +i z -) 



-iAz sign (3?Az) + D(k + k') 2 



(6) 



We have denoted Az = z + — z-, np the density of 
states at the Fermi energy, and D the diffusion con- 
stant. This form of the low-energy behavior holds also 
for stronger scatterings on impurities. Multiple scatter- 
ings only renormalize characteristic parameters A and 
D of the singular part of the electron-hole vertex for 
Az — > 0.0 Equation © is a generalization of the dif- 
fusion pole from real to complex frequency differences. 

To demonstrate incompatibility of the Ward iden- 
tity with causality of the self-energy we define a real- 
frequency function 

AW(u>) = I £[E*(k, E — id) — £*(k, E + u)]. (7a) 



From the Ward identity we obtain its representation via 



the irreducible electron-hole vertex 



kk' 



-Af$(E - u, E)] %G R (k' : E) + A R £(E + u, E) 
x [G R [E + u)-G R , (E)] - A R £ (E-u,E) 

x [G2(E- U )-G%(E)]} . (7b) 

We have used superscripts R, A to denote the complex 
half-planes of the energy arguments (upper, lower) from 
which the real value is approached. Analyticity of the 
self-energy as a function of energy E inside the energy 
bands induces analyticity of the function AW(u>), i. e., 
all its derivatives are continuous functions and bounded 
at u> = 0. 

The vertex function A RA is not, however, analytic and 
its nonanalyticity is transferred to AW(u>). It is straight- 
forward to evaluate AW smg (u>), the leading singular part 
of AW(u>) in the limit u> — » 0. We rescale the momenta 
q — > x = qyj D/\lj\ in A RA and perform momentum in- 
tegration in a d-dimensional space. The leading singular 
behavior is determined by the small-momentum asymp- 
totics of the first term on the r.h.s. of Eq. (|7b(l with the 
singular electron- hole vertex from Eq. © for 77 = + . We 
explicitly obtain 
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for d ^ 41, 
for d = 4/, 



(8) 



where K and kp are a dimensionless constant of order 
one and the Fermi momentum, respectively. |lf| Equa- 
tion (|7af) . however, dictates that AWT ng {u)) = almost 
everywhere inside the energy bands. We hence have a 
contradiction to the premise of validity of the Ward iden- 
tity Eq. (J2J for real frequency differences. Q. E. D. 

We can see from Eq. iJSJ that the nonanalyticity in 
AW(lo) strongly depends on the spatial dimension d. It 
has a square root character in odd dimensions, contains a 
sign function in dimensions d = 2(21 + 1) and a logarithm 
in d = 41. The nonanalytic behavior is smoothened with 
increasing dimensionality of the system, i. e., the higher 
the dimension the higher is the derivative of AVF(w) dis- 
playing an (infinite) jump at u> — 0. It is the [(d— 2)/2]th 
derivative, where [x] is the smallest integer greater than 
or equal to x. As expected, there is no nonanalyticity in 
the mean-field, d — 00, limit of the electron-hole vertex 
and one cannot use a 1/d expansion around the CPA in 
two-particle functions. 

The nonanalyticity in AW(u>) demonstrates that we 
cannot simultaneously keep causality and the Ward iden- 
tity. Causality of the self-energy (reality of eigenenergies 
of the random Hamiltonian) cannot be broken. We hence 
are forced to admit a violation of the Ward identity for 
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real frequency differences z + — z_ = 2uj =/= in theo- 
ries with a causal self-energy and the Cooper pole in the 
electron-hole irreducible vertex. 

Giving up the Ward identity has quantitative con- 
sequences for the diffusive behavior of disordered elec- 
trons. Most importantly, the asymptotic form of the dif- 
fusion/Cooper pole deviates from Eq. JfjJ. This devia- 
tion in the low-energy asymptotics of the electron-hole 
correlation function can be explicitly manifested in high 
spatial dimensions. There we expect only a weak, negli- 
gible momentum dependence of the self-energy as well as 
of integrals of the vertex function. The Cooper singular- 
ity in the vertex function is sufficiently smoothened by 
momentum integration. If we introduce a parameter 
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E fl (k,B + U )-i^A^(£ + W) B) 



(G R (k' 7 E + cj)-G A (k',E)) 



(9) 



j=0 



we can represent the leading low-energy asymptotics of 
the zero-temperature electron-hole correlation function 
in high spatial dimensions as 



$* A 



^0B + W ,B;q) 
2irnp 



N 2 
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The parameter A increases with the disorder strength. 
In particular, it is the imaginary part of the frequency 
derivative of the electron- hole vertex, 3<9A e/l /du>, that 
drives the system toward the Anderson localization tran- 
sition, where it eventually diverges. The weight of the dif- 
fusion pole expressing the number of electron 
states available for diffusion, decreases when approaching 
the localization transition as well as an effective diffusion 
constant D/(l + A). They both vanish at the Anderson 
localization transition. There is hence no diffusion pole 
in the localized phase. 

Violation of the Ward identity can be consistently ex- 
plained and does not actually contradict any fundamen- 
tal physical law. The Ward identity holds only within 
the Hilbert space spanned by extended Bloch waves, the 
eigenstates of the nonrandom kinetic-energy operator. [j| 
Violation of the Ward identity then expresses the fact 
that we cannot construct a single Hilbert space for all 
configurations of the random potential. There are new, 
disorder-induced and configuration-dependent localized 
states orthogonal in the thermodynamic limit to Bloch 
(plane) waves. The statistical weight of these states is 
not negligible and hence the state vector of the electron 
in random media has a non-zero projection onto a space 
of localized states at any disorder strength. That is why 
the number of electrons in the space of extended states is 
not conserved when the disorder strength increases. The 



weight of the localized states in the decomposition of the 
electron state vector increases and reaches its saturation 
at the Anderson localization transition, where the state 
vector becomes orthogonal to Bloch waves. 

To conclude, we proved in this Letter that the Cooper 
pole, Eq. (JSJ, in the electron-hole vertex leads in sys- 
tems with disordered electrons to deviations from the 
Ward identity. The nonanalytic behavior due to the 
Cooper pole in the electron-hole vertex is not sufficiently 
smeared out by momentum integration in the Ward iden- 
tity to produce a self-energy analytic in its energy ar- 
gument. The observed deviation from the Ward iden- 
tity in causal theories is explained by nonconservation of 
electron states available for diffusion, i.e., plane waves. 
Due to the existence of disorder-induced localized states, 
the electrons may escape from the Hilbert space of ex- 
tended states into the space of localized ones, orthogonal 
to Bloch waves. 
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